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Abstract 

We study quantum mechanical and classical stability properties of Reissner- Nordstrom 
deSitter (RNdS) spacetimes, which describe black holes with mass M and charge Q in a 
background with cosmological constant A > 0. There are two sources of particle production 
in these spacetimes; the black hole horizon and the cosmological horizon. A scattering 
calculation is done to compute the Hawking radiation in these spacetimes. We find that the 
flux from the black hole horizon equals the flux from the cosmological horizon, if and only 
if \Q\ = M, indicating that this is a state of thermodynamic equilibrium. The spectrum, 
however, is not thermal. We also show that spacetimes containing a number of charge 
equal to mass black holes with A > 0, have supercovariantly constant spinors, suggesting 
that they may be minimum energy states in a positive energy construction. As a first 
step in this direction, we present a positive energy construction for asymptotically deSitter 
spacetimes with vanishing charge. Because the construction depends only on a spatial 
slice, our result also holds for spacetimes which are asymptotically Robertson- Walker. 
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1. Introduction 



In this paper we will discuss quantum mechanical and classical aspects of the stability 
of Reisner-Nordstrom deSitter (RNdS) spacetimes, in which the magnitude of the electric 
charge Q is equal to the mass M. In the remainder of the introduction, we give the 
argument, based on euclidean quantum field theory, that \Q\ = M RNdS black holes 
are stable end points of the process of Hawking evaporation. In section 2, we present a 
lorentzian scattering calculation of particle production in RNdS spacetimes, which shows 
that the flux of particles across the cosmological horizon is equal to the flux across the black 
hole horizon, if and only if Q = M. Such a calculation is useful for a number of reasons. 
As well as confirming the Euclidean argument, in this case, it provides information on the 
spectrum of particles produced. Moreover, generically a spacetime will not have a real 
euclidean section, and we will be limited to lorentzian techniques. The multi-black hole 
generalizations of the \Q\ — M RNdS spacetimes |I| are examples of nonstatic spacetimes 
which should have interesting thermodynamic properties^. In section 3, we show that 
\Q\ = M RNdS spacetimes admit spinor fields which are constant with respect to a certain 
super-covariant derivative operator. The choice of derivative operator is motivated by 
supergravity and should lead to a version of the positive energy theorem relevant to RNdS 
black holes 0. Here, we present a more limited, preliminary result relevant to positivity of 
a suitably defined energy in asymptotically Robertson- Walker spacetimes without charge.^ 

RNdS spacetimes have metric and gauge field given in static coordinates by 



where the cosmological constant A is assumed to be positive. For a range of values of 
Q and M, the spacetime has three Killing horizons; inner and outer black hole horizons 
and a cosmological horizon. The extremal limit, in which the inner and outer black hole 
horizons become coincident, occurs for M < \Q\, with equality in the case A = 0. 

There are two sources of particle production in a RNdS spacetime; the black hole 
horizon and the deSitter horizon || . It is interesting to ask whether these two sources can 
ever be in a state of thermal equilibrium. The deSitter horizon would then be like the wall 
of a box containing a black hole in thermal equilibrium with a bath of radiation. Several 
facts suggest that such an equilibrium exists when the magnitude of the charge \Q\ and 
mass M of the black hole are equal and furthermore that the equilibrium is stable, as 
we now review. 

One can assign temperatures Tbh and Ta s to the black hole and deSitter horizons 
based on the periodicity that must be imposed on the imaginary time coordinate in order 

3 The euclidean picture of black hole thermodynamics can also be misleading, as in the case 
of charged dilaton black holes || 

4 It has come to our attention that work on a positive energy construction in asymptotically 
deSitter spacetimes has recently been done by T. Shiromizu. 



ds 2 



f{R)dT 2 + f(R)- 1 dR 2 + R 2 dfl 2 1 A 
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to obtain a regular euclidean section of the metric ([[]) at the given horizon [ffl. The 
temperatures defined in this way are given by 

2irkT bh = K bhl 2nkT ds = K dSl (2) 

where k is Boltzmann's constant and n b h and KdS are the magnitudes of the surface gravities 
at the black hole and deSitter horizons respectively. For \Q\ = M, one finds K b h = K dS 
and the two temperatures are equal. Further, if M > \Q\ then one finds that K b h > KdS 
(and vice versa), suggesting that the black hole radiates (or absorbs) energy to reach 
equilibrium at \Q\ = M. This is also the configuration of maximal gravitational entropy. 
The gravitational entropy is given in terms of the areas A bbi and A ds of the black hole and 
deSitter horizons by S = 4 (A b h + A^s)- For an infinitesimal perturbation between RNdS 
solutions with fixed charge the first law of thermodynamics [|Jstates that — KdS^^-dS = 
K bh 5A bh . i Hence 5A tot = 5A bh (l - ^) = for n bh = n dS . 

These arguments based on euclidean quantum field theory suggest that if quantum 
mechanical processes are taken into account, so that for example the area of the black hole 
horizon can decrease, then a charged black hole with \Q\ ^ M will evolve to a state with 
\Q\ = M. In this paper, we give the results of a Lorentzian scattering calculation of the 
particle production in Hawking radiation. We compare the flux of particles crossing the 
cosmo logical horizon at late times ( i.e., particles "coming from the black hole") to the flux 
which crosses the black hole horizon (i.e., particles "coming from the deSitter horizon"). 
We find that the fluxes are equal if and only if \Q\ = M. In this sense, a \Q\ = M black hole 
is in equilibrium with the deSitter background, in agreement with the Euclidean picture. 

2. Particle Production 

We compute the rate of particle production in the RNdS spacetimes (1), for a massless 
scalar field cf>, by finding the mode mixing coefficients and Ptjfui' 0- ^ ne conformal 

diagram for the relevant portion of RNdS is shown in fig. [I]. The region is bounded by 
the white hole, black hole, past and future deSitter horizons, which have replaced past 
and future null infinity in the case with zero cosmological constant. The first step is 
to define what is meant by particles, i.e. we must choose a time coordinate near each 
boundary, which defines the positive frequency modes there. Define the following Kruskal 
type coordinates which are well behaved near the appropriate horizons 

1 „ .. „ 1 

(3) 



U hh = e~ KbhU , Vhh = e KbhV 



'bh — e , Vbh 

K>bh K>bh 

U dS = — e^ su , V dS = -—e~^ v 



KdS K d S 

where, in standard notation u and v are defined in terms of the tortoise coordinate r* 
according to 

dv 

u = t - r* ,v = t + r* ,dr* = —. (4) 

f 2 



5 This can be generalized to include arbitrary, not necessarily stationary, stress energy pertur- 
bations, by the methods of ||. Then there appears an additional volume integral of the perturbed 
energy density on the right hand side. 
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Near the black hole horizon, the metric then has the limitting form 



ds 2 ps K bh dU bh dV hh , (5) 

showing that these are good cooordinates near the black hole horizon. One can similarly 
show that the coordinates are also good near the deSitter horizon. 

For A = 0, one defines positive frequency near I~ by </> ~ e~ luJV , and positive frequency 
on the white hole horizon by <fi ~ e~ lujJJbh . This choice of positive frequency on the white 
hole horizon correctly gives the outgoing thermal flux at I + , allowing one to replace the 
collapsing body in Hawking's original calculation with this boundary condition on the 
field (see e.g. This suggests that for RNdS spacetimes, positive frequency at each 

horizon should be defined with respect to the appropriate Kruskal coordinate there. One 
can check that this is equivalent to choosing the affine parameter on e.g., the white hole 
horizon to define positive frequency there. 

The Klein-Gordon equation for (p near any of the horizons reduces to the free wave 
equation. As in the case with no cosmological constant, the potential term due to the 
background gravitational field decays exponentially in r* near a horizon. Consider then 
a pure positive frequency outgoing wave near the deSitter horizon at late conformal time, 
^ ~ e~ iujUdS (see fig. |l]). In the geometries optics limit (as with A = 0), finding the form 
of this wave propagated back to the white hole horizon reduces to finding the dependence 
of the coordinate Uds on the coordinate Ubh- Using ||), it then follows that on the white 
hole horizon 

,e-^ 2 ^\ U bh <0 
GUU bh )=(f>UUds(U bh ))~ { (6) 

0, U bh > 

where rj = Kds/ K bh and £ 2 = '^(^^) r> - The particle production follows from determining 
the negative frequency portion of this wave on the white hole horizon, 



bh 

ui.uj' 



0J 



J dU bh e- luJ ' Ubh GUU b h) 



(7) 



Similiarly, there is emission "from" the deSitter horizon as seen by an observer out- 
side the black hole horizon at late times. Consider a positive frequency wave which is 
entering the black hole horizon (see fig. I]), 0^ ~ e ~ lujV bh_ Then in the geometries optics 
approximation, on the past deSitter horizon, the wave is given by 

FUVds) = MV bh (V dS ) ~ { 6 " S ' dS< (8) 

0, V dS > 

where fi 2 = 1/ K b hi)-I K ds) * ■ Similarly to ([F[) , fl^f u , is given in terms of the fourier transform 
of (H). For general values of Q and M, the functions and G w appearing in (H) and @ 
are related according to 

G u (x) =F^(x" 2 ). (9) 
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We see that the two functions are equal for r\ = 1, which occurs when \Q\ = M . We then 
have the main result of this section 



ft!,% if and only if \Q\ = M. (10) 

This implied that for each horizon, the flux of particles absorbed is equal to the flux of 
particles absorbed. 

The spectrum of emitted particles is given by = f (Zu/|/3 W)a /| 2 . We can estimate 
the above integrals using the stationary phase approximation. It is simpler to work with 
the coefficients a W;LU ' = —0uj,-u>'- The stationary phase approximation can be used to 
evaluate this integral, and (3 WjU i' is then gotten by analytically continuation, after noticing 
that implies that ol w ,<jj' is analytic in the lower half ui' plane. 

For the case \Q\ = M, when the surface gravities are equal, one finds in the stationary 
phase approximation 

\3 J 2 = ^ 1 e ->*T /a Mi) 

This approximation can be checked, because in the \Q\ = M case the integral can be done 
exactly. One can show that a u ^i is only a function of the combination x = ^l^w'l 1 / 2 . 
The function a(x) can then be shown to satisfy the differential equation 

a"(x) + -a' + (4 - -^a = 0. (12) 

x x 2 

Hence a is proportional to the Bessel function Ji(x). Again analytically continuing, we 
then have /5 Wia> / ^ |wu/| 1//2 ). Using the asymptotic behavior of the modified Bessel 

function Ki(x), this agrees with our previous expression for large values of the argument 
x 1 and fixes the constant of proportionality to be V2/k. For small values of the argument, 
one then finds |/3 w ,w'| 2 ~ 2 ' • However, the geomtrics optics approximation which was 
used in doing the calculation means that the results are only valid at high frequencies. We 
would expect that the wavelength of the propagating wave must be small compared to the 
length scales and r^s- Using a lower limit of 2-k /r^ in the integration over frequencies, 
one finds for the spectrum 

Nu w —e ^7S (13) 
4a; 

This is not a thermal spectrum, so the situation is probably best described as a thermo- 
dynmic equilibrium. If the lower cutoff were taken to be proportional to u, rather than 
l/ r Wi, the spectrum would be thermal. 

There are several limits one can take in order to check this result. Letting k — ► 
above, corresponds to keeping \Q\ = M and letting the cosmo logical constant A approach 
zero, so that the spacetime approaches extremal Reissner-Nordstrom. In this limit goes 
to zero, as it should. Secondly, one can consider setting Q = 0, and then letting A — >• 0, 
so that the metric approaches Schwarzschild. This is equivalent to letting r^s approach 
zero, while keeping finite. The particle production (^) from the black hole can again 
be evaluated in the stationary phase approximation. One finds that the coefficients a w>U} ' 
approach those for a Schwarzschild black hole. 
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3. Supercovariantly Constant Spinors for A > 

We now turn to issues relevant to the classical stability of \Q\ = M RNdS spacetimes. 
First, recall that for A = a generalization of the positive energy theorem [fT0|, ||ll| implies 



that the ADM mass M of an asymptotically flat spacetime satisfies the inequality 

M>\Q\, (14) 

where Q is the total electric charge of the spacetime. This inequality is saturated for 
spacetimes which admit a spinor field e satisfying V M e = 0, where V M is a certain super- 
covariant derivative operator arising naturally in ungauged N = 2 supergravity. In the 
context of supergravity V^e gives the variation of the gravitino field under an infinites- 
imal local supersymmetry transformation. Spacetimes having supercovariantly constant 
spinors have unbroken supersymmetries. Amongst asymptotically flat electrovac solutions, 



the bound ([14]) is saturated by the Majumdar-Papapetrou multi-black hole solutions ||12|| , 
which represent collections of extremal (|Q| = M) charged black holes. These spacetimes 
each admit two supercovariantly constant spinors. Calculations of the effective potential 
in supergravity have shown that these spacetimes remain ground states at the quantum 
level as well Jl]| . 



A nonzero cosmological constant A = —3g 2 arises in gauged N = 2 supergravity, where 
g is the coupling constant of the gravitino with the U(l) gauge field. If the coupling g is 
real, as one would require in a quantum theory, then the cosmological constant is negative. 



Romans [|14| has written down the supercovariant derivative operator V M , which arises 
in gauged N = 2 supergravity, with real coupling g and showed that \Q\ = M RNadS 
spacetimes each admit two supercovariantly constant spinors@. At the classical level one 
can also consider imaginary values of g and hence positive values of the cosmological 
constant. With the substitution g = iH, the supercovariant derivative operator is given 
by 

V M e = f 4 V p + HA» + % -H llM + % -F pa Yl a l^ e (15) 

We are then interested whether or not \Q\ = M RNdS black holes admit spinor fields 
satisfying V^e = 0. If this is the case, then it seems promising that a version of the 
positive energy theorem can be proved, showing the classical stability of these spacetimes. 

In fact we can find supercovariantly constant spinors for a more general class of space- 
times, which includes the \Q\ = M RNdS spacetimes. These are the deSitter analogues of 
the MP multi-black hole solutions, which were found in PJ. In spatially flat Robertson- 
Walker type coordinates, these spacetimes are given by 



ds 2 = - -pJ^ro dt 2 + a(t) 2 ^l(x) 2 5 ii dx i dx j , A = ^-dt, 

O = 1 + V , ^\ , , a(t) = exp(Ht). 
~ a\x — Xi 



(16) 



6 These spacetimes have naked singularities, because for A < the extremal limit occurs with 
\Q\ < M. 
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We find that these spacetimes each have two supercovariantly constant spinors, which are 
given by 

I p -Ht/2 

e=-=e = ^^e (17) 
v afi v O 



where e is a constant spinor satisfying the projectionB. (1 — irf)e = 

4. Positive Energy Constructions 

As stated above, the supercovariant derivative operator ([15]) should be the starting 
point for a positive energy construction, showing that a suitably defined mass for asymp- 
totically deSitter spacetimes is minimized for fixed charge by the spacetimes (|T6|). The 
relevant mass in the asymptotically deSitter case should turn out to be that given by Ab- 



bott and Deser in [15]. In the present note, we will address only the case Q = 0. Such 



a construction for A < is given in | T1| . We will return to the Q ^ case in a future 
publication [Q. The result we present here, however, will be valid for spacetimes which 
are asymptotic to spatially flat Robertson- Walker, not necessarily asymptotically deSitter. 
This is because the construction takes place only on a spatial slice, and a spatial slice of 
Robertson- Walker can be embedded in a deSitter spacetime. In this case the asymptotic 
mass will be related to the boundary term given in |J, in the context of peturbations of 
Robertson- Walker spacetimes, and discussed further in (16[ . 

Let the metric g^ v of a spacetime asymptotically approach the metric for a Robertson- 
Walker spacetime g^ on each t = constant spatial slice, where 

ds 2 RW = g™dx^dx u = -dt 2 + a 2 {t)5 l3 dx i dx j . (18) 
With vanishing gauge field the supercovariant derivative operator (|15]) becomes 



V At e = V M e+~#(£) 7a e, (19) 

where H{t) = a{t)/a{t) may depend on time. 

When the metric is exactly Robert son- Walker, there exist four supercovariantly con- 
stant spinors e G satisfying V M e D = 0. This can be seen as follows. In five dimensional 
Minkowski spacetime there are four constant spinor solutions to V M "i/> = 0. Projecting 
this equation onto the four dimensional deSitter hyperboloid, one finds that there are four 
solutions to the equation V^V = (V M + ^lal^)^ = 0. Letting e = (1 — ij^)ip, we havei 
V^eo = 0. A t = constant spatial slice in a general Robertson- Walker spacetime ( [18| ) 
can be embedded in the deSitter hyperboloid with the appropriate value of H{t). On a 
given spatial slice then, there are four solutions to Vje = , where j is a spatial index. 
Explicity, these are given by 

e G = (1 - z 7 S )(l - | * W + 7 S ))e (20) 



7 In our conventions hatted indices denote frame indices and unhatted indices are coordinate 
indices. Hence {t a ,7 P } = 2^ and {7", 7"} = 2g^ '. 

8 Note that (j 5 ) 2 = +1, so this is not a projection. 
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where e is any constant spinor. 

Since the positive energy construction will take place on a spatial surface, this is 
actually sufficient to construct a positive energy statement, for any spatial slice in the 
coordinates of equation (|18|). However, there are solutions which satsify the full equation 
V M i? = 0. These are conveniently given by E = exp(± J ^dt')e , where 7*60 = =Re D . 

The latter projection is equivalent to (7* ± 7 s ) e = 0. 

Let Qij and be the metric and extrinsic curvature on a spatial slice, which approach 
the metric and extrinsic curvature of a Robertson- Walker spatial slice with Hubble constant 
H. Assume that gij and Kij solve the Einstein constraint equations with sources and 
2|£. Let e be a solution to the contracted equation 

y V 3 -e = (21) 

on this spatial slice ( where the index j runs over spatial values only) with the boundary 
condition that e approaches one of the constant spinors e a at large r. Let V be a volume 



in the spatial slice. Then following the analysis of JL7fl we find the following relation 

/ da^V'e = [ d 3 x y/g U^{T U - Tg w + l£ 7 V)e + (V^V'e) , (22) 
Jdv Jv v 7 

where T~ w = 3H(t) 2 /8ir is the energy density of a Robertson- Walker spatial slice. If 

black holes are present in V, then the boundary dV has components at the black hole 

horizons. These boundary terms can be made to vanish by imposing a suitable projection 

on e, as shown in ||11|| . This leaves the outer boundary of V, which we will take to be 



spatial infinity. In the case A = 0, this would be related to the ADM 4-momentum (T7j . 
The volume integral on the right hand side of (|22|) can, in the context of asymptotically 
Robertson- Walker spacetimes, have both positive and negative contributions. The second 
term, involving derivatives of the spinor field, is positive definite. However, the first term, 
which is a difference between the matter stress energy and the background Robertson- 
Walker energy density, can be either positive or negative, depending on whether there are 
mass overdensities or underdensities. If the matter stress energy satisfies the dominant 
energy condition, however, this term cannot be more negative than minus the energy 
density of the background itself. Furthermore, for asymptotically deSitter spacetimes it is 
natural to consider this difference to be positive. 

One expects that the boundary term at spatial infinity gives an analogue of the ADM 
4-momentum for asymptotically Robertson- Walker spacetimes. Such a quantity has been 



defined for asymptotically deSitter spacetimes by Abbott and Deser |15[ . For each asymp- 
totic killing vector in an asymptotically deSitter spacetime, one can construct a conserved 
charge given by a two dimensional boundary integral on a spatial slice. These charges 



all vanish for pure deSitter. Hence, as is argued in the anti-deSitter case |pj| , we expect 
that for asymptotically deSitter spacetimes, the boundary term in (|22] ) gives the Abbott- 
Deser mass. In the example of Schwarzchild-deSitter, one can evaluate the boundary term 
explicity, and it gives the mass parameter. 



More generally, the McVittie metric |18| describes a mass embedded in a Robertson- 
Walker spacetime, 

1 m_ m 

ds 2 = (- 2 -^) 2 dt 2 + a 2 (t)(l + —)% j dx t dx j . (23) 

1 + 7^ 2ar 

Zar 
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In this case, the exact solution to the Dirac equation (21) can be found. The boundary 
term can then be evaluated and shown to give the mass parameter m. Hence, if a metric 
approaches (^3|) at large radius up to terms of order 1/r 2 , the boundary term will continue 
to give the appropriate mass parameter. In principle, one can find a general expression for 
the boundary term in ( |22|) by solving (|2l| ) asymptotically on a spatial surface to sufficiently 
high order, as in [17]. In practice, here, this is complicated by the nonzero background 
extrinsic curvature and we have not yet found a general expression. We expect that the 
final result will be related to the boundary integral which arises in the analysis of || on 
perturbations of Robertson- Walker spacetimes. 
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Figure Captions 



Fig. 1. Shown is the conformal diagram for the section of RNdS relevant to the particle 
production calculation. The region is bounded by the white hole, black hole, past 
deSitter, and future deSitter horizons. Also indicated are the choices for positive 
frequency waves at each horizon. To compute p^^i , one studies wave propagation 
along a null geodesic Uds = constant, parallel to the blackhole horizon. The wave 
is taken to be positive frequency and outgoing as it crosses the future deSitter 
horizon. Propagating back to the white hole horizon, the wave is decomposed 
into positive and negative frequency parts with respect to the coordinate Ubh 
on the white hole horizon. Near the white hole horizon the wave vanishes for 
Ubh > 0, since that is inside the black hole. Similiary, to compute (3^,, one 
considers a wave propagating on a null geodesic Vbh = constant, parallel to the 
future deSitter horizon, for a wave which is positive frequency and ingoing at the 
black hole horizon. 
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